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Assuming the standard cosmological model as correct, the average linear size of galaxies 
with the same luminosity is six times smaller at z = 3.2 than at z = 0, and their average 
angular size for a given luminosity is approximately proportional to z" 1 . Neither the 
hypothesis that galaxies which formed earlier have much higher densities nor their lumi- 
nosity evolution, mergers ratio, or massive outflows due to a quasar feedback mechanism 
are enough to justify such a strong size evolution. Also, at high redshift, the intrinsic 
ultraviolet surface brightness would be prohibitively high with this evolution, and the 
velocity dispersion much higher than observed. We explore here another possibility to 
overcome this problem by considering different cosmological scenarios that might make 
the observed angular sizes compatible with a weaker evolution. 

One of the models explored, a very simple phcnomenological extrapolation of the 
linear Hubble law in a Euclidean static universe, fits the angular size vs. redshift depen- 
dence quite well, which is also approximately proportional to z~ 1 with this cosmological 
model. There are no free parameters derived ad hoc, although the error bars allow a 
slight size/luminosity evolution. The type la supernovae Hubble diagram can also be 
explained in terms of this model with no ad hoc fitted parameter. 

WARNING: I do not argue here that the true Universe is static. My intention is 
just to discuss which theoretical models provide a better fit to the data of observational 
cosmology. 

Keywords: Cosmology: miscellaneous; galaxies: statistics. 



1. Introduction 

The analysis of the dependence of the angular size of some sources with redshift 
was for many decades one of the most important geometric tests of cosmological 
models. Different cosmologies predict different dependences for a given linear size 
and this can be compared with the data from observations. The test, first conceived 
by HoyleP, is simple in principle but its application is not so simple because of the 
difficulty in finding a standard rod, a type of object with no evolution in linear size 
over the lifetime of the Universe. 



1 
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It is well known that the application of the angular size (9) vs. redshift (z) 
test gives a rough dep endence of 9 oc z~ x for QSOs and radio gal axies at r adio 
wavelengthJ^EEEGOZf f or f[ r st ranked cluster galaxies in the Optica j8 | 9 | 10 | ll | 

and 

for the separation of brightest galaxies in clusters or in QSO-galaxy pairs of 
the same redshiftP^. The deficit of large objects at high redshifts with respect to 
the predictions of an expanding Universe is believed to be an evolutionary effect 
by which galaxies were smaller in the past (e.g.,^j), or a selection effect (e.g., ^^J). 
Thus, the 9 oc z _1 relationship, as predicted by a static Euclidean Universe, would 
be just a fortuitous coincidence of the superposition of the 9{z) dependence in the 
expanding Universe and evolutionary/selection effects. 

Som e other studie s have tried to find better standard rods. Ultra-compact radio 
sources ^^^^^^ were used to carry out an angular size test: a dependence 
different from 9 oc z" 1 and closer to the predictions of expanding Universe models 
was found. The test was even used to ascertain not only whether or not the Universe 
is expanding but also to constrain the different cosmological parameters. However, 
these applications are not free from selection effect and, as will be discussed in 
£15. 2[ interpretation of the results of these tests is not so straightforward. 

Another proposaPS used the rotation speed of high redshift galaxies as a stan- 
dard size indicator since there is a correlation be twe en size and rotational velocity 
of galactic disks. This method was indeed appliecPH over a sample of emission-line 
galaxies with 0.2 < z < 1, with the result that Einstein-de Sitter cosmology is 
excluded to within 2-cr, and that small galaxies (with fixed rotation velocity lower 
than 100 km/s) should show a strong evolution in size (at z = 1 a factor two smaller 
than at z = 0) and no evolution in luminosity within the concordance cosmology, 
while large galaxies do not evolve significantly, either in size or in luminosity. Un- 
fortunately, this method requires spectroscopy, so the sample has an upper limit of 
z ~ 1 even with very large telescopes, and the uncertainties are so large that they 
cannot be interpreted directly without certain assumptions conc erni ng evolution 
models. In principle, looking at figures 4 and 5 of Marinoni et aiPU, one sees no 
reason to exclude a 9 oc z _1 dependence. 

The aim of this paper is to repeat the angular size test for galaxies within a wide 
range of redshifts (z = 0.2 — 3.2) in optical-near infrared surveys (equivalent to 
the optical at rest). Data from high spatial resolution surveys available nowadays, 
such as those carried out with Hubble Telescope or FIRES, provide useful input 
for this old test of th e angular size with new analyses and interpretations. Recent 
analysis of these 

dat£ f22l23l22 

has shown that the linear size of the galaxies with 
the same luminosity should be much lower than locally. However, this is true only 
if we considered the standard cosmological model as correct. In the present paper, 
I will do a reanalysis of these data and consider different cosmological scenarios, 
which will shed further light on the degeneracy between expansion + evolution and 
non-expansion. 
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2. Data 

Angular effective radii, defined as circular ized S ersic half-light radii within which 

50% of the light is present, are taken from 1221231 f or galaxies with 0.2 < z < 1 and 
1041 

from Trujillo et al¥^tov galaxies with z > 1. Both samples separate approximately 
early-type and late-type galaxies by means of the exponent (ns) of their Sersic 
profile: ns > 2.5 for early-type galaxies and ns < 2.5 for late type galaxies. 

| 22 | 23 | provide data and angular size measurement of the GEMS survey with 
two Hubble Space Telescope colors (F606W and F850LP). In total, they have 929 
galaxies. The data, processing and photometry are discussed by Rix et a/PS. 

Trujillo et erfJ" use near-infrared FIRES data of z > 1 galaxies in the HDF-S 
and MS 1054-03 fields to derive the angular size in the rest-frame V filter. In total, 
there are 248 galaxies with 1 < z < 3.2. There are 14 more galaxies with z > 3.2, but 
they are very few, very luminous, and their photometric redshift determination was 
not very accurate; indeed, the available on-line data through the FIRES Website, 
http: / /www.strw.leidenuniv.nl/~fires| gives new recalculated redshifts and some of 
them are very different. The data processing and photometry are discussed in detail 
by Labbe et alW for HDF-S and Forster-Schreiber et alW& for the MS 1054-03 
field. 

From these galaxies, I take only those with 3.4 x 10 10 < Lv,rest(Lo.v) < 2.5 x 
10 11 , a total of 393 galaxies (271 galaxies from GEMS with z < 1, and 122 galaxies 
from HDF-S/MS 1054-03 fields with z > 1). The lower limit is the same as that 
adopted by Trujillo et ai.^l avoiding the faintest ones in order to have a more 
homogeneous sample. The maximum limit is to avoid the galaxies away from the 
range of local galaxies for which the relationship between radius and luminosity was 
explored (explained in <J3]). This is an almost complete sample for redshifts up to 
« 2.5 [24] an( j i s i ncom plete for 2.5 < z < 3.2. As we will see later, it is unimportant 
whether the sample is complete or not since the test is independent of the luminosity 
of the galaxies, but within this restriction I will concentrate on the analysis of the 
brightest galaxies. A higher limit at the lowest luminosity (the sample would be 
complete up to z = 3.2 for Ly, rest > 6.7 x 10 10 L Q ) would reduce the number of 
galaxies too severely and the statistics would be poorer. In any case, I will also 
comment on the results for these higher luminosity lower limits (see §3.2|) . 



3. Angular size test 

I am going to analyze the variation of the angular size, 9y,rest, of the galaxy rest- 
frame V with the redshift. The angular size in the z < 1 sample was measured 
at A = 9450 A (filter F850LP), 6\ , which, due to the color gr adient s, is slightly 
different from the angular size at V-rest. This difference is smafPSESl DU t I apply 
the following correction to it: 
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'V.rest 



1 _ o.ii ( *vO±»2 _ i 
1 _ 0.08 ( Ay( , 1+z) - 1 

\ Ao 




(1) 



with Ao = 9450 A, and Ay = 5500 A. In any case, this correction is very small (less 
than 4%) and the results would not change significantly if it were not applied. In 
HDF-S/MS 1054-03 data, the size was already measured in the filter which gives V 
at restt^H so no correction is necessary. 

Since we have a wide range of luminosities and types of galaxies, there is a huge 
dispersion of sizes for a given redshift, with Malmquist bias, but this dispersion and 
bias can be reduced by defining 9* as the equivalent angular size if the galaxy were 
early-type with a given V-rest luminosity (I take 10 10 L Qi y; however, the variation 
of this number does not affect any result, just the calibration in size): 

_ Q R(10 10 L Q ,y,ng>2.5) 
°* = Wf 7T\ ' ^ ' 

where R(Lv.rest,ns) is the average radius of a galaxy for a given luminosity and 
exponent (ns) of the Sersic profile. The number ns is affecte d b y an important 
uncertainty for galaxies with very small angular size (6 < 0.125" ) , which produces 
some extra dispersion. In any case, the dispersion of R values is moderate^, and, 
given that 6*(z) does not contain the dispersion of luminosities, it will present a 
much lower dispersion than 9(z). Certainly, 9*(z) contains the spread of 9 and R so 
the dispersion due to random errors in these quantities is larger for 9* than for 9; but 
the dispersion due to the spread of luminosities dominates, so, as said, will present 
a dispersion much lower than 9{z). This definition also avoids selection effects due 
to Malmquist bias, since 9* for a given redshift should be nearly independent of the 
luminosity of the galaxy, at least on average. 

Shen et aZ. 28 (Eqs. 14-15) give the median radius of a galaxy of given r' SDSS - 
luminosity (K-correction applied) and ns for local SDSS galaxies. In total, the radius 
(measured for z w 0.1 in the r-band, which is more or less equivalent to V-band at 
rest) as a function of the absolute magnitude in r' SDSS -iest M r >_sDSS is: 

R(M r ,-sDSS,n s ) (3) 

10 -0.260M r ,_ SDSS -5.0 6) ns > 2.5 \ 

10 -0.104M r ,_ SDSS -1.71[ 1 + 10 -0.4(M r ,_ SDSS +20.91)]0.25 ; ns < 2 .5 J kpC 

We are not considering the evolution in this Eq. ([3]); all discussion of the effects 
of evolution will be considered in 21 To translate this relationship into a V-rest 
luminosit y, w e must make a color correction, as in Mcintosh et al¥^ [however Tru- 
jillo et alV^ interpolate between the rest- frame g-band and r-band]. Taking into 
account that ((V — r' SDSS AB )) =0.33 for early-type galaxies and = 0.30 for late- 
type galaxieJ^ (this already includes the transformation of Vega to AB system; no 



2, 2010 16:14 WSPC/INSTRUCTION FILE sizes 



Angular size test 5 

evolution is considered here), and that My.Q = +4.79 (Vega system) we will have 
(Ly in units of 10 10 L 0! y) 



DfT . f A L° V 65 , n s >2.5\ , 

" s) = i B L™(1 + fliLv) - 35 , n s < 2.5 } kpC ' (4) 

with Ao — 1.91, -Bo = 2.63, Si = 0.692. These parameters are valid assuming the 
concordance cosmological model (Hq — 70 km/s/Mpc, Q, m — 0.3, ft\ = 0.7). For 
other cosmologies, we have to calibrate the relationship between luminosities and 
radii with the corresponding luminosity and angular distances for zsdss ~ 0.1. See 
Table Q] for the values of Ag, Bq and B\ with other cosmologies. This relationship 
of Shen et is fitted with galaxies of —19 > M r > —24 for early types and 

-16 > M r > -24 for late types, which is equivalent to 2.5 x 10 9 < L v < 2.5 x 10 11 
L Q y for early types and 1.6 x 10 8 < Ly < 2.5 x 10 11 L© y for late types. Our 
galaxies are within these limits. 

In order to derive the luminosity, Ly trest , for a given redshift (z) and the rest 
flux for the filter V, Fy trest , we need the distance luminosity di,{z) from different 
cosmologies, such that (without considering neither evolution nor extinction) 



V 4irFy rest 

And the average equivalent angular size evolution with z should be compared 
with the prediction of the same cosmology, which, in principle, without taking into 
account the evolution, should be given by: 

dj^^z) ^* (6) 

Q*,pred{z*) 

where dA{z) is the angular distance, and i?* is the equivalent physical radius of the 
galaxy associated with the equivalent angular size that is, if the galaxy were 
early-type of V-rest luminosity 10 10 L Qj y. 



3.1. Different cosmological scenarios 

(1) Concordance model with Hubble constant Hq = 70 km/s/Mpc, Vl m = 0.3, 
fl A = 0.7: 



d A ( z ) = f dx , (7) 

M > H (l + z)J ^(I + ^ + Oa' K> 



d L {z) = {l + z) 2 d A (z). 



(8) 
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(2) Einstein-de Sitter model [Eq. © with ft A = 0, Q m = 1]: 

1 



d A {z) = 



H (l + z) 



1 - 



VT+z. 



(9) 




d L (z) = (l + z) 2 d A (z). (10) 

Although this is not the standard model nowadays, there are some researchers 
who stil l consider it more appropriate than the concordance model (e.g., 
| 30 | 31 | 32 fy. About the compatibility with type la supernovae data, see discus- 
sion in ^5.31 

Since the distances for objects at z = 0.1 are 0.943 times the distances of the 
concordance model, the corrected relationship between luminosity and radius 
is Eq. with A = 1.90, B = 2.61, B 1 = 0.699. 
(3) Friedmann model of negative curvature with f2 = 0.3, Q\ = 0, which implies a 
term of curvature fl K = 0.7 [23, Eq. (26)]. 



d A (z)= TT - ~~r ^ sinhl y/Sl K I m = d Q X _ = ) , (11) 



d L {z) = (l + z) 2 d A (z). (12) 

I will use this model of an open universe to check that a model different to a 
flat universe does not change the results significantly. Since the distances for 
objects at z = 0.1 are 0.970 times the distances of the concordance model, the 
relationship between luminosity and radius with corrected calibration is Eq. (}4|) 
with Ao = 1.93, B = 2.59, B x = 0.736. 
(4) Quasi-Steady State Cosmology (QSSC), Q m = 1.27, fl A = -0.09, O c = -0.18 
(C-field density) EH 

d A (z) = —4- r / =, (13) 

H (i + z)Jo ^/n c (i + x) 4 + n m (i + x) 3 + n A 

d L (z) = {l + z) 2 d A {z). (14) 

This cosmology is n ot the standard model, but it can al so fit many data on 
angular size test J^l or Hubble diagrams 

for SNe I^5|36|37| 

The expansion with 

an oscillatory term gives a dependence of the luminosity and angular distance 
similar to the standard model, adding the effect of matter creation (C-field) with 
slight changes depending on the parameters. The parameters of this cosmology 
are not as well constrained as those in the standard model. He re, I use the 
best fit for a flat (K = 0) cosmology given by Banerjee et oZ.^I Cl m = 1.27, 
= —0.09, Vl c = —0.18, which corresponds to r] = 0.887 (amplitude of the 
oscillation relative to 1), xq = 0.797 (ratio between actual size of the Universe 
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and the average size in the present oscillation) and maximum allowed redshift 
of a galaxy z max = 6.05 (note however that the ma xim um observed redshift has 
risen above 8 nowadays according to some authors, ^Sl). Other preferred sets of 
parameters give results that are close. T he values most used are K = 0, ft\ — 
-0.36, T) = 0.811, z max = 5 |35 | 36 | 37 | 39 | which imply Qm = x 63 ^ c _ _ .27, 

but I avoid them because they do not allow galaxies to be fitted with z > 5. 
Parameters with a curvature differe nt fr om zero (K ^ 0) also give results that 
are very close in the angular size tesiP^l 

Since the distances for objects at z = 0.1 are 0.943 times the distances of 
the concordance model, the relationship between luminosity and radius with 
corrected calibration is Eq. (gj with A = 1.94, B = 2.56, B x = 0.778. 
(5) Static euclidean model with linear Hubble law for all redshifts: 



d A {z) = -jjrZ, (15) 
B<3 



d L { Z ) = y/T+zd A (z). (16) 

These simple relations indicate that redshift is always proportional to angular 
distance, a Hubble law. We assume in this scenario that the Universe is static; 
the factor y/1 + z in the luminosity distance stems from the loss of energy 
due to redshift without expansion. There is no time dilation, and precisely be- 
cause of that the factor is not (1 + z). The caveat is to explain the mechanism 
different from the expansion/Doppler effect, which gives rise to the redshift. 
This cosmological model is not a solution which has been explored theoreti- 
cally/mathematically. However, from a phenomenological point of view, we can 
consider this relationship between distance and redshift as an ad hoc extrapola- 
tion from the observed dependence on the low redshift Universe. Our goal here 
is to see how well it fits the data, and forget for the moment the theoretical 
derivation of this law. 

Since the luminosity distance for objects at z — 0.1 is 0.966 times the luminosity 
distance of the concordance model and the angular distance is 1.115 times the 
angular distance of the concordance model, we must adopt approximately Eq. 
flU with A = 2.23, B = 2.98, B 1 = 0.741. 
(6) Tired- light/simple static euclidean model : 



d A (z) = -^-Hl + z), (17) 



d L (z) = VT+ld A (z). (18) 



This is again a possible ad hoc phenomenological representation which stems 
from considering that the photons lose energy along their paths due to some 
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interaction, and the relative loss of energy is proportional to the length of that 
path (e.g. 021), i.e. 

d 4- = ——E. (19) 
dr c 

Of course, as in the previous case, this ansatz is very far from being considered 
as the correct one by most cosmologists, but it is interesting to analyze its 
compatibility with the angular size test too. 

For the calibration of Eq. (£[]), I use the fact that the luminosity distance for 
objects at z = 0.1 is 0.921 times the luminosity distance of the concordance 
model, and the angular distance is 1.063 times the angular distance of the 
concordance model, so A Q = 2.26, B = 2.92, B\ — 0.816. 
(7) Tired-light/ "Plasma redshift" static euclidean model: 



d A (z) = ^-ln(l + z), (20) 

d L (z) = (l + z)V 2 d A (z). (21) 

The plasma redshift application^^ §5. 8) used dL{z) = (1 + z) 3 ^ 2 dA(z) instead 
of dh{z) = (1 + z) 1 / 2 d y i(z) to take into account an extra Compton scattering 
which is double that of the plasma redshift absorption. For the calibration of 
Eq. Q: A, = 2.00, B = 2.78, B x = 0.674. 

The volume element in a static uiverse is different from the volume element in 
the standard concordance model. Particularly for the standard concordance model 
the comoving volume element in a solid angle du and redshift interval dz is (for null 
curvature, which is the case of the concordance model) 



dVr 



concordance 



. 3 

C 



r 

Jo 



tlx 



v /n m {i+x) 3 +n A 



while for the two first static models it is 



-dujdz (22) 



c 



C^static-lin.Hub. = [ Jf ) ^dudz, (23) 



^static-tir.1. = ( [hl ^ + f 2 du>dz. (24) 

\H Q J (l + z) 

Hence, if we wanted to evaluate the evolution of some quantity per unit comoving 
volume for the static universes, we must multiply the result in the concordance 
model by the factor rfv " 1 ?,"- cordi " lcc . 

J "^static 
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3.2. Results 

The results of the test are plotted in Fig. [1] for the concordance model, with the 
equivalent angular size of each individual galaxy, and Fig.[2]for the different models, 
with the representation of the average value of log 10 in bins of A log 10 (z) = 0.10. 
I do a weighted linear fit in the log-log plot. Since I am calculating an average of 
the logarithm for the angular sizes, the possible error in individual galaxies should 
not significantly affect the value of the average. The error bars in each bin of the 
plot are the statistical errors. Note that there are values of 0* lower than 0.03", but 
they do indeed correspond to measured values of 9 > 0.03"; 9* is lower than 9 at 
the highest redshifts because the luminosity in V-filter of those galaxies is higher 
than 10 10 L Qi y. 

As can be observed, the average equivalent angular size gives a good fit to 
a 0*(z) — Kz~ a law, with values of K, a given in Table [2] No fit is totally in 
agreement with the cosmological prediction without evolution and extinction. The 
static model with a linear Hubble law is not very far from being compatible with the 
data: I get 9*(z) = 0.136z -097 , near the expected dependence (9*(z) — 0.109z _1 ) 
although with a slightly larger size. 

NOTE: the normalization of the angular size in any model prediction ( solid 
line in Fig. 0) stems from the Shen et all ^ calibration with the corresponding pa- 
rameters Aq, Bq and B\ in each cosmology. It does not stem from a fit. 

If we separate elliptical (ns > 2.5) and disk galaxies (ns < 2.5) for the con- 
cordance model, the plots in Fig. [3] are obtained, with a higher slope for elliptical 
galaxies (a = 1.18 ± 0.04) than for disk galaxies (a = 0.80 ± 0.07). This might be 
due to a different evolution of disk and elliptical galaxies, but they are likely to be 
due to different systematic errors for elliptical galaxies and disk galaxies due to sys- 
tematic errors in the value of n$- The galaxies with 2.0 < ng < 2.5 are suspected of 
being strongly contaminated by ellipticals since the obtained ns tends to be lower 
than the real one ({231, Fig. 1). If we take disk galaxies only with ng < 2.0, then 
a = 0.92±0.06, closer to unity. On the other hand, elliptical galaxies with 9 < 0.125" 
are also strongly contaminated by disk galaxies (J^ 1 , Fig. 2) which are compacted 
by a wrong measure of ns and consequently give a smaller radius than the real one. 
If we take elliptical (ns > 2.5) galaxies only with 9 > 0.125": a = 1.03 ± 0.05. 
Therefore, from the present analysis and within the systematic errors, we cannot be 
sure that the angular size test gives different results for elliptical and disk galaxies. 
However, when all the elliptical and disk galaxies are put together, the excesses and 
deficits more or less compensate; there are approximately as many elliptical galaxies 
misclassified as disk galaxies as disk galaxies misclassified as elliptical galaxies. See 
further discussion on the systematic errors in §3.31 

In Fig. 01 I analyze the dependence of a on the luminosity of the galaxies, and 
we see that there is no significant dependence, a = 1 is more or less compatible 
with all the subsamples of different luminosity. 
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Cosmology: concord, std. £2 =0.3, £1 =0.7 



0.1 : 



* Data (GEMS) 

□ Data (HDF-S; MS 1054-03) 



1001 



025 



0.5 



I 



redshift (z) 



Fig. 1. Log— log plot of the equivalent angular size (0*) vs. redshift (z) in the concordance cos- 
mology. Left: (z < 1) GEMS data. Right: (z > 1) HDF-S and MS 1054-03 data. The average of 
log 10 6>* in bins of A log 10 (z) = 0.10 is represented with asterisks and squares with statistical error 
bars. 



Values of the parameters Aq, Bq and B\ in the relationship between radius and 
luminosity of a galaxy, Eq. calibrated with different cosmological models. 



Cosmology 


A 


Bo 


B x 


Concord. fi m = 0.3, il\ = 0.7 


1.91 


2.63 


0.692 


Einstein-de Sitter 


1.94 


2.56 


0.778 


Friedmann f2 TO = 0.3 


1.93 


2.59 


0.736 


qssc n m = i.27, n A = -o.o9, n c = o.is 


1.94 


2.56 


0.778 


Static linear Hubble law 


2.23 


2.98 


0.741 


Static, tired light/simple 


2.26 


2.92 


0.816 


Static, tired light/plasma 


2.00 


2.78 


0.674 


St. lin. Hub. law, ext. a v = 1.6 x 10~ 4 Mpc" 1 
St. tired light, ext. a v = 3.4 x 10" 4 Mpc -1 


2.21 
2.22 


3.03 
3.01 


0.696 
0.719 
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Weighted fit of the average equivalent angular size for the data in Fig. [2] to a law 
6*(z) = Kz~ a . Error bars give only statistical errors and do not include systematic 
errors as explained in ^3-31 



Cosmology 


K 


a 


Concord. TO = 0.3, il\ = 0.7 


0.1251 ±0.0041 


0.957 ±0.045 


Einstein-de Sitter 


0.1721 ±0.0054 


0.834 ±0.044 


Friedmann VL m = 0.3 


0.1412 ±0.0046 


0.956 ±0.045 


qssc n m = i.27, n A = -0.09, n c = o.is 


0.1810 ±0.0057 


0.817 ±0.044 


Static linear Hubble law 


0.1363 ±0.0044 


0.969 ±0.045 


Static, tired light/simple 


0.2132 ±0.0066 


0.717 ±0.043 


Static, tired light/plasma 


0.0915 ±0.0031 


1.177 ±0.047 


Concord., early type (ns > 2.5) 


0.1059 ± 0.0044 


1.181 ±0.043 


Concord., late type (ns < 2.5) 


0.1441 ±0.0070 


0.805 ±0.067 


Concord., early type (ns > 2.5), 6 > 0.125" 


0.1226 ±0.0060 


1.031 ±0.051 


Concord., late type (ns < 2.0) 


0.1600 ±0.0073 


0.915 ±0.057 


Concord., L v > 6.8 x 10 10 h Qy 


0.1315 ±0.0083 


0.995 ±0.098 


Concord., L v > 1.02 x 10 11 L QyV 


0.1106 ±0.0117 


1.313 ±0.186 



St. lin. Hub. law, ext. a v = 1.6 x 10" 4 Mpc" 1 0.1055 ± 0.0035 1.113 ± 0.046 
St. tired light, ext. a v = 3.4 x 10~ 4 Mpc" 1 0.1565 ± 0.0050 0.801 ± 0.044 



3.3. Selection effects, and errors in the angular size measurement 
and luminosity-radius relationship 

As said previously, the definition of 6* avoids selection effects due to Malmquist 
bias, since 6* for a given redshift should be nearly independent on how luminous 
the galaxy is, at least on average. 

Nonetheless, although i?„ = 0*cU is independent of the luminosity and type of 
the galaxy at z = 0, its possible evolution might depend on both parameters and 
this would give different average 0* when the selection of galaxies is different. In 
any case, whatever the sample of galaxies is, if we are going to interpret the factor 
between the cosmological prediction and the observed average 6* as a product of size 
evolution, this factor would reflect the average evolution of i?» in that sample. In 
our case, we are observing the factors for the average population with 3.4 x 10 10 < 
Lv,rest(L(?)) < 2.5 x 10 11 . For the concordance model the value of i?» at z — 3.2 is on 
average ~ 6 lower than i?* at z — 0. At low redshift we know that the relationship 
of Eq. dU is correct (hence, the average size i?* should not depend on the luminosity 
of these galaxies). Therefore, the ratio of sizes between high/low redshift objects 
will depend only on the variation of Eq. (j4} in galaxies selected at high redshift. We 
can say that the average size of the selected galaxies at z — 3.2 is w 6 times lower 
than the galaxies at low redshift with the equivalent luminosity and Sersic profile. 
In other words, if there is a (very strong) evolution in the size of the galaxies for a 
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Fig. 2. 

Log-log plot of the average of log 10 0*, where 9„ is the equivalent angular size, vs. 
redshift (z). Bins of Alog 10 (z) = 0.10. Error bars only represent statistical errors; 
for the systematic errors, see text in §3.31 The seven plots are for the seven different 
cosmologies described in §3.11 Solid lines are the model predictions (the normaliza- 

inoi 

tion stems from the Shen et al¥ l21 calibration with the corresponding parameters 
A , B and B\ in each cosmology). Dashed lines are the best weighted linear fits. 
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Fig. 3. Log-log plot of the average of log 10 where 0* is the equivalent angular size, vs. redshift 
(z). Bins of A log 10 (z) = 0.10. The two plots are for the concordance cosmology separating elliptical 
galaxies (rig > 2.5) from disk galaxies (ng < 2.5). 




Fig. 4. Log— log plot of the average of log 10 0* , where 0* is the equivalent angular size, vs. redshift 
(z). Bins of Alog 10 (z) = 0.20 and Alog 10 (z) = 0.30 respectively. The two plots are for the 
concordance cosmology only for galaxies with Ly > 6.8 X 10 10 Lq y (122 galaxies) and Ly > 
1.02 X 10 11 L Qj v r (48 galaxies) respectively. 



fixed luminosity (or a variation in luminosity for a fixed radius), the factor by which 
the galaxies are smaller will depend on which galaxies are selected, but in any case 
this factor will represent an average galaxy shrinking factor. An average factor of 6 
in size reduction will mean that there are galaxies with size reduction by a factor 
larger than 6, an d o ther galaxies with size reduction by a factor smaller than 6. 

Trujillo et a/!2^J(§4.3) discussed the robustness of the "average" angular size 
measurement, whether it is affected by some other biases or selection effects. Their 
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tests showed that the results, as presented here, are more or less robust and the 
difference in the results if we apply some minor corrections due to uncertainties or 
incompleteness in the radius distribution of the sources are A((log 10 R)) ~ 0.1, thus 
A0*/0* ~ 20% . 

There is a minimum angular size that Hubble telescope or FIRES can observe 
below which the uncertainties of 8 and ns are very high, with important systematic 
deviations. Normally, the angular sizes are overestimated for these low angular sizes, 
so the observed average 0* increases with z as more and more of the smaller angular 
size galaxies are included in the average. Moreover, the very compact galaxies are 
not included in the sample because they might be classified as stars. This affects 
mainly the z > 1 data, where the limit for low systematic errors is around = 0.125" 
1231. Suppression of these galaxies would change the average of 9* at z > 1 by 10- 
30%. Systematic errors for 6 < 0.125" are up to 50% or lower ( 24 , Fig. 2), which 
produces a systematic error on the average 9* lower than « 15%. 

Another effect to investigate is the error in the relationship of Eq. (j4j. This 
represents the average radius for a given luminosity and there is some dispersion 

IOQI 

of values with respect to it (given by Shen et al¥^). I am not interested here in 
these statistical errors of R because this would just affect a dispersion of values of 
9* without changing its average value. Our concern is about the systematic errors 
in Eq. (j4}. I will analyze five sources of systematic errors: 

(1) The color correction of ((V— r' SDSS AB ))=0.33 for early-type galaxies and = 0.30 
for late-type galaxies might change if there is a preferential galactic type 
within the group of early- or late-type galaxies is given for some redshift. Indeed, 
Fukugita et al^^ give a value of ((V — r' SDSS AB )) of 0.36 for E type, 0.31 for 
SO, 0.32 for S a b, 0.29 for S& c , 0.28 for S c d- That is, there are variations up to 
« 0.03 with respect to the average in early types, and 0.02 for late types. In the 
worst cases, these systematic errors in colors would produce a systematic error 
of m 2% in 6*, which is negligible. 

(2) A misclassification of an elliptical galaxy as disk galaxy or vice versa produces 
an error up to a factor 2 in the value of R derived from Eq. (j4|), which leads 
to an error of factor < 2 in 9*. The effect of the uncertainty in ng for small 
angular sizes has already been checked by Trujillo et aZ.2%4.3) and the error 
of the average of 9* is lower than ~ 20% for the sample with all of the galaxies. 
However, when I divide the galaxies into elliptical and disk types (see [3]) the 
effect might be larger, as was noted previously concerning Fig. |3l 

(3) Systematic errors in the lumin osit y at V-rest due to errors in the photometric 
redshift amount less than 4% SB, which implies errors in 9* less than 2.5%. 
However, the systematic error of the luminosity measurement in faint elliptical 
galaxies amounts around 15% 123 which translates into 6% of error in 9*. 

(4) The calibration of eq. (|4| (calibration of the parameters Aq, Bq and Bi) is done 
for cosmologies different to concordance model assuming that the SDSS galaxies 
have redshift zsdss ~ 0.1. There is indeed a dispersion of redshift values in 



2, 2010 16:14 WSPC/INSTRUCTION FILE sizes 



Angular size test 15 



SDSS galaxies around this average. If we set zsdss = 0.15 we would obtain, for 
instance for the static linear Hubble law ansatz, Aq = 2.32, Bq — 3.15, B\ — 
0.721; the variations in the average measured are negligible (<< 1%), but 
the predictions of the cosmological model for 0* (proportional to A$) increase 
by w 6%. That is, the ratio between measured and predicted decreases by 
6%. An error of 0.05 in the average redshift is among the worst cases, so we can 
say that the systematic error in the calibration should be less than « 5%. 
(5) Another systematic error comes from the V m ax correction of the selection effects 
in the SDSS data made by Shen et a/P^(§3.1). Shen et aL^lmake a correction 
of selection effects by assigning a weight to each galaxy inversely proportional 
to the maximum comoving volume within which galaxies identical to one under 
consideration can be observed^. This correction mainly affects faint elliptical 
galaxies, as can be observed in fig. 2 of Shen et aZpS and amounts to <^ 
0.7 x exp(— 0.7Ly) in the average radius of the galaxies for elliptical galaxies of 
luminosity Ly (in units of 10 10 L S) y). 

The weights depend on the cosmological model, since the comoving volumes at 
a given redshift change in the different cosmological models. Moreover, there is 
an intrinsic systematic error in using the standard model: Shen et a/PSl used 
the integration of physical volumes instead of comoving volumes [see Eq. (8) of 
Shenei a/.^j, which leads to systematic errors of V max up to ~ 50% [there is a 
factor (1 + z) 3 between comoving and physical volumes]. In order to calculate 
the effect of a change of cosmology, apart from the recalibration of Ao, Bq 
and Bi, I would need to have all their SDSS data at hand and repeat the full 
analysis of their fit with the new conditi ona l maximum volume V max values for 
each cosmology [see Eq. (9) of Shen et al^\ using the correct comoving volume 
Vmax instead of the integration of physical volumes, which is not possible for 
us. Assuming a total systematic error in the volume of AV max ~ 0.5V^* ax , 



Since Ly > 3.4 in our selected sample, this may justify systematic errors in 
R(L) up to 3% . 

Summing up, apart from the statistical errors plotted in Fig. [2J there are sys- 
tematic errors in the average that may amount up to 30-40% for the general case 
with all the galaxies 

3.4. Relationship with the surface brightness test 

The average surface brightness of a galaxy with total flux F\, rest and half-light 
circularized angular size 9 is: 



A((log 10 R)) <~ 0.15 x exp(-0.7{L v )), 



(25) 



SB = 




(26) 



7T(9 2 
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Using eqs. ([5]), ([6]); and = (l+z) l ^ 2 dA, with i — 1 if static, and i = 4 if expanding, 
it is 

The intrinsic surface brightness [without the (1 + z) % dimming factor] is 

SB ° = W- (28 » 

Given the relationship between radius and luminosities of eq. Q, we find that the 
average intrinsic surface brightness in V-rest (assuming no size evolution) should 
follow: 



-0.30 



-M-/1/.1 ■/>-) = { ^' L o y 4 8 " 5 > 2 ' 5 > • (2!)) 

87r 2 Bg(l+SiL v ) - 50 ' < 2.5 

The surface brightness for a given luminosity should be independent of the rcdshift 
if there is no evolution in size. I can also define an equivalent surface brightness to 
avoid the dependence on the luminosity: 

SB (1Q W L QiV ,> 2.5) 10 10 L QiV . . 

*~ 27T02 SB (L v , rest ,n s ) ~8^Rl(l + z y [6U > 

This is indeed the surface brightness test, or Tolman testP^l. In Fig. [2] for the 
static models, it is observed that the data of 0*(z) without size evolution are fitted 
more or less (within the statistical+systematic error). Thus, i?» is nearly constant 
in a static model for all redshifts and we obtain (SB*) oc (1 + z)^ 1 , while for the 
expanding model we would need a strong evolution in R*(z) and the average surface 
brightness would be (SB*) oc i?(z)~ 2 (l + z)~ 4 . In order to give the same (SB*)(z) 
as in the static Universe, assuming that d^z) is approximately the same (which is 
nearly true for the concordance and the linear Hubble law luminosity distances; see 
3H1J), R(z) w R(z = OKI + z)" 3/2 . 

Lubin & Sandage^ obtained (SB) oc (1 + z) 1 with i = 1.6 — 3.2 for a sample 
with z < 0.9 depending on the filter (R or I), and the initial hypothesis (expanding 
or static). Lerner351 obtained i = 1.03 ± 0.15 for z <^ 5, with data in wavelengths 
from ultraviolet to visible from Hubble Space Telescope (see also LernerS^l), while 
Nabokov & Baryshe\E2l obtained i = 3 — 4 with the same type of data but without 
including K-corrections. Andrews^l obtained i = 0.99 ± 0.38, i — 1.15 ± 0.34, with 
two different samples of cD galaxies. Lubin & Sandage argue that their data are 
compatible with an expanding Universe if we take evolution into account. Lerner 
criticizes Lubin & Sandage for not using the same range of wavelengths at rest for 
the high and low redshift galaxies but instead using K-corrections with many free 
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parameters, which is less direct and more susceptible to errors. Moreover, Lerner 
argue against the evolution that the intrinsic ultraviolet surface brightness of high 
redshift galaxies would be extremely large, with impossible values (see E )5.4.ip . 

4. Evolution of galaxies in expansion models 

From the plot in Fig. concord, model), we see that, in order to make the concor- 
dance model compatible with the data on angular size, we must assume an evolu- 
tion such that the galaxies at high redshift are much smaller than at low redshift 
(z <^ 0.2). For instance, at redshift z = 2.5, the galaxies should be on average 
« 4.6 times smaller. Trujillo et dP^ obtained a lower average factor (~ 3) because 
they used only galaxies with the constraint 8 > 0.125"; however in the error bars of 
Trujillo et al. result, this bias effect of avoiding 8 < 0.125" galaxies is included, and 
within this error bar it is compatible with our result. The m ost massive galaxies 
are thought to be contracted by a factor of ~ 4 up to z = 1.5 S3 and of 5.5 up to 
z = 2.3^, but Mancini et alA^ think that the extra compactness of these galaxies 
can be understood in terms of fluctuations due to noise preventing the recovery of 
the extended low surface brightness halos in the light profile, so those factors might 
be not real. For galaxies with redshift z = 3.2 in our Fig. [2j concord, model), the 
average ratio of sizes would increase to a factor 6.1 (comparing the linear fit with 
the theoretical prediction). Separating by types, the ratio would be 4.5-6 for disk 
galaxies and 6-10 for elliptical galaxies, ta king into account th e systematic uncer- 
tainties commented in £ 33.21 Other author 

j5() | 51 | 52 | 53 | 54 | 56 | 3 31 

find similar results 

at high redshift. 

If we have galaxies that are on average 6 times smaller than the equivalent 
galaxies (of same type and luminosity) at low redshift, this means that the V-rest 
luminosity density (inversely proportional to the cube of the size) of these galaxies is 
R3 200 times higher than at low redshift. The surface brightness in V-rest is increased 
by a factor ~ 40. Is this situation possible? 

4.1. Effect of the expansion of the Universe 

In the standard picture of the expanding models, the expansion does not affect 
the galaxies because there is no local effect on particle dynamics from the global 
expansion of the universe: the tendency to separate is a kin ema tic initial condition, 
and once this is removed, all memory of the expansion is loslP-H Note, however, that 
there are other views on the topic of whether there is expansion of gal axies due to 
the expansion of the Universe and the nature of the expansion itself EHES] rj ar k 
energy or cosmological vacuum might hav e so me influence on the size of the galaxies 
BUM!! but this effect would be small. Lee^ suggested instead that the size of the 
galaxies increases as the scale factor of the universe assuming dark matter models 
based on a Bose-Einstein condensate or scalar field of ultra-light scalar particle, 
which is another heterodox idea to explore. Nonetheless, apart from this kind of 
proposals, within standard scenarios, the galaxies do not expand with the Universe. 
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An indirect way in which the expansion has an effect is in the formation of 
galaxies. In the theoretical ACDM hierarchical scenarios, galaxies formed at higher 



redshift should be denser^ since, to decouple from expansion, structures must have 
a given density ratio with t he surrounding density, which is larger at higher redshift. 
Some authors (e.g., |24 | 50 |j nave use d this argument to explain the apparent size 
evolution. However, the observed redshift in galaxies is not its formation redshift, 
so the application of this idea is not straightforward. As a matter of fact, the stellar 



populations of most local massive elliptical galaxies are very oloP^, and formed 
before the age corresponding to z ~ 3, so we cannot say that galaxies now are 
larger because of this effect since most of them were formed > 10 Gyr ago. The 
difference of formation epoch of the galaxies observed now and at z ~ 3 is not 
large. Moreover, the hierarchical scenario in which massive galaxies form first do 



not represent the observed Universe appropriately, as I argue in 

Furthermore, the theoretical claim by Mo et aZp^J derived from the models that 
galaxies which formed earlier are denser is not in general observed. If it were true, 
we should observe that at low redshift the youngest galaxies (formed later) should 
be much larger for a given mass than the oldest galaxies of age 12-13 Gyr. There is 
already evidence that this is not the case: the densest galaxies are young instead of 



okPl And we can check with our own sample within 0.2 < z < 3.2 that the color 
of elliptical galaxies is not correlated with size: Fig. [5] Redder elliptical galaxies are 
older and, for a given redshift, indicate earlier formation, which should be equivalent 
to smaller size, at least statistically. This correlation is not observed at all: linear 
fits in the four redshift ranges of Fig. [5] all give slopes compatible with zero within 
la except for the range —0.4 < log 10 (z) < —0.1, which gives 'jnf^rS~ = 
+1.0 ± 0.5, a 2cr correlation but in a direction opposite to prediction that galaxies 
are larger when redder=older (formed earlier). 

Therefore, it is not a question of comparing the formation of galaxies at dif- 
ferent redshifts but the evolution of galaxies already formed, either isolated or in 
interaction/merging with other systems. 



4.2. Younger population of high redshift galaxies 

The main argument in favour of the evolution in size for a fixed luminosity is that 
the younger a galaxy is the brighter it is, and we expect to see younger galaxies 
at high redshift. Therefore, galaxies with radius smaller than i?„ in the past will 
produce the same luminosity as galaxies with that radius at present. How much 
brighter? 

Using Vazdekis et aZ.'iP^ synthesis model, with revised Kroupa IMF, we can 
derive the mass-luminosity ratio in a passively evolving elliptical galaxy as a func- 
tion of its intrinsic (B-V) color and its luminosity. The metallicity degeneracy is 
broken with a n it erative method which uses the relationship between stellar mass 
and metallicity^^. With this method, the average mass-luminosity ratio of ellipti- 
cal galaxies (not affected by extinction) at the last bin of Fig. |2] {z = 2.5 — 3.2) is 



2, 2010 16:14 WSPC/INSTRUCTION FILE sizes 



Angular size test 19 



q>2.5 







1 


□ □ 

□ D n □ 








« „ !" », 

' °° "Ink - f ' 


1 














, J $1.°"* V ° " ' 






■t 


■ • * "".i ■ * , 


< 






* * * " " 


« 0.1 






* * 






o -0.7<log 10 z<-0.4 








o -0.4<log 10 z<-0.1 




0.01 




, -0.1<log 10 z<0.2 








, 0.2<log |0 z<0.5 








1 






0.2 0.4 0.6 0.8 



Fig- 5. 

Equivalent linear size (normalized with Shen et al.'s^^ calibration) vs. (B — V) color 
at rest for elliptical galaxies in the concordance cosmology. 

(M*/L v ) = 0.5 [N = 6, (L v ) = 10 11 L , ((B - V)) = 0]. Rudnick et aV^{%AA) 
analysis of the same galaxies also obtain a quite similar mass-luminosity ratio. For 
the elliptical galaxies with (Ly) > 5 x 10 10 L Q of the two first bins (z = 0.20 — 0.32): 
(M*/L v ) = 2.2 [N = 7, (L v ) = 7 x 10 10 L , ({B - V)) = 0.94], that is, a mass- 
luminosity ratio 4.4 times larger. Assuming a variation of the luminosity linearly 
dependent on the time, the variation of this ratio would be by a factor 5.8 between 
z=0.1 (the average redshift of SDSS galaxies, which are the reference of the size cal- 
ibration in Shen et a/PHJ) and z = 3.2. This is an acceptable estimation for elliptical 
galaxies. Kauffmann et a ;(S8l(Fig. 13) showed that SDSS late-type galaxies have a 
m ass- luminosity ratio around 2 times lower than early-type ones, and Rudnick et 
aZ3^(Fig. 9) showed that blue galaxies at z sa 3 also have mass-luminosity ratios 
around 0.2-0.3, so we also keep this number of ~ 6 for disk galaxies. This factor of 6 
might be affected by important errors (see §5.4p . and could be much lower, though 
not much larger. 

Given that we are comparing galaxies with the same luminosity, the galaxies at 
z = 3.2 would have 6 t ime s lower stellar mass than at z = 0.1. Hence, from Eqs. 
(17)-(18) of Shen et a/PSl, we find that galaxies at z = 3.2 should be 2.7 times 
smaller than at z = 0.1 if they are elliptical, or 2.0 if they are disk galaxies. These 
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factors are much lower than the measured values of 6-10 and 4.5-6 respectively. 
A factor in size 2-4 for elliptical or 2-3 for disk galaxies, including this luminosity 
e volu tion correction remains (in rough agreement with the results by Trujillo et 
aZ.^j). Therefore, the argument of "younger population in higher redshift" does not 
serve to justify the present data. In order to explain the observed size increase in 
terms only of luminosity evolution of the stellar population, we would need to set 
M^/Ly with respect to the SDSS galaxies a factor 25-60 for elliptical galaxies, or 
50-100 for disk galaxies, too much! 



4.3. Mergers 

If the luminosity density is not due to an increase of the luminosity of each star, it 
must be due to a redistribution of the mass density to make it more compact. Why? 
Expla nations in terms of mergers pr oliferat e in the literature (see discussion by Refs. 
I ^4 | 47 | anc j references therein). Refs. I ^4 | 47 | SU pp 0r t the merger scenario calculating 
the size vs. mass ratio and observing how it also decreases with redshiflQ, but this 
calculation also depends on the cosmological model used and can give different 
results for different cosmologies. Furthermore, Khochfar et alP^ point out that the 
presence of higher amounts of cold gas at high redshift mergers of ellipticals also 
produces size e volu tion. 

Trujillo et al\^i cites the paper of Boylan-Kolchin et aZp^as a possible powerful 
mechanism to increas e th e radius in massive elliptical galaxies. Each merger would 
follow a law R oc M^EQIwith a = 1.0 for a pericenter of 15 kpc and lower for higher 
pericenter distances or higher for lower pericenter distance. So when two galaxies 
of the same luminosity approach each other reaching a minimum distance of 15 
kpc and merge, the tot al r adius will be ~ 2.00 times the radius of each individual 
galaxy, while Shen et a/PSlfor local galaxies give a — 0.56: a factor « 1.47 in radius 
on average when we double the mass. This means that each major merger (fusion 
of galaxies of the same luminosity) gives an extra factor of 1.36 for ellipticals in 
angular size with respect to eq. (j4]) relationship. We would need an average of 2.3- 
3.6 (= 36) ' i- e - 1-36 2-3-3 ' 6 = 2 — 3) major mergers along the life of each elliptical 
galaxy to justify a factor 2-3 in radius. 

The observed rate of mergers i s in deed much lower than the necessary rate 
to justify these numbers. Lin et a/P^ with statistics of close galaxy pairs (r < 
20/i~ 1 kpc) up to z=1.2 show that only ~ 9% of the luminous galaxies (—21 < 
Mb < — 19) would have a major merger during their lives since z = 1.2. De Propis 
et al^^ get similar merger ratios for local galaxies (z < 0.12) than Lin et alU^ 



a In Refs.ElMH the masses were estimated from the colors of the galaxies assuming sola r m etallicitv 
for the determination of the mass— light ratio and a Salpeter IMF model. Trujillo et aXfi^ multiply 
the masses by a factor 0.5 to correct for the difference of K roupa and Salpeter IMF, which is only a 
very rough approximation. The masses of Trujilloei alW^ are more accurately calculated, with an 
uncertainty of a factor two. All these assumptions introduce significant errors into the calculation 
of the mass-light ratio. Therefore, we must take these results with care. 
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for 0.5 < z < 1. The number of cumulative mergers increases up to 22% if we 
allow larger pair separation (r < 30/i _1 kpc) and up to 54% ESI if we allow a broader 
definition of m ajor merger including pairs of galaxies with mass ratio up to four. 
Ryan et aZp^ calculated a number of 42% galaxies undergoing some merger up to 
z = 1 for more luminous galaxies (Mb > —20.5) [and 62% for all z], also for mass 
ratios up to four, and r < 20h~ 1 kpc. In our case, comparison with major mergers 
of equal mass galaxies and with average pericenter around 15 kpc should be made, 
so a number of 10-20% (up to z = 3.2 is ~ 50% larger than up to z — 1) of mergers 
would be the amount to compare with the 2-4 mergers (200-400%) we need for 
the size evolution. Or we may account for mergers with mass ratios up to 4 and 
pericenter radii larger than 15 kpc on average, although these would not produce an 
increase of an extra factor of 1.36 for ellipticals in angular size with respect to the 
Eq. (j4j relationship but lower. Therefore, we would have to compare the number 
of ~ 50% with the large number of mergers of this kind to produce the observed 
evolution in radius (~ 10 mergers per galaxy; probability of merging ~ 1000%). 

Moreover, I suspect that the number of mergers is overestimated. First because 
not all galaxy pairs become mergers. Second, because of the contamination of in- 
terlopers. There is an uncertainty of radial distance due to proper motion of the 
galaxies apart from the Hubble flow, and many of the identified pairs of galaxies in 
the projected sky are not real pairs in 3D space. There may be a chanc e su perposi- 
tion in the line of sight of galaxies. The situation is worst for Ryan et a/f^l, who use 
spectrophotometric redshifts; and interlopers do not introduce a statistical error, as 
they say, but a systematic one. However, the order of magnitude of the pairs with or 
without interloper s sh ould be more or less the same, the interlopers being at z < 3 
lower than ~ 30% Also, De Ravel et alU^, for instance, using spectroscopically 
confirmed pairs, get similar merger ratios. Third, because timescales of mergers 
increa se s lightly with redshift and are longer than assumed in most observational 
studies^. 

About the CAS method of estimation of merger rat ios ba sed on the identifi- 
cation of major mergers with highly asymmetric galaxiesPEi two things may be 
commented: 1) The authors consider that all asymmetries in principle associated 
with starbursts are due to major mergers, but the mechanisms which trigger impor- 
tant amounts of star formation might be different from major (ratio 1:1) mergers; 
in particular, minor mergers may also trigger asymmetric star formation. Strongly 
disturbed systems, indicative of recent strong interactions and mergers, account for 
only a small fraction of the total star formation rate density^! 2) Interlopers, galax- 
ies and stars with very different redshifts projected as background or foreground 

objects in the line of sight produce an important amount of a ppa rent distortion in 

iVqi 

the galaxies, especially at high redshift, where Conselice et aly-^ claim that a high 
fraction of galaxies are major mergers. We must bear in mind that Hubble images 
may detect very faint galaxies and there are more than 5 million galaxies per square 
degree with m z < 30 pAl extrapolated from the counts up to magnitude 28), one 
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galaxy in each square of 1.6" x 1 .6" on average, which, mixed with the main galaxy 
(m z < 27 in Conselice et aiJ^j), produce apparently distorted galaxies. Also, De 
Propis et a/.^l checked that the contamination is happening in low redshift galaxies 
with foreground stars. Moreover, we see the galaxies once they are "presumably" 
merged but we do not see enough galaxy pairs at 1 < z < 2 when both galaxies 
of the mer ger are separated (e.g., according to the already overestimated ratio by 
Ryan et alJ^ 1 ). This indicates at least that the asymmetries are not produced by 
major mergers but by minor mergers (ratio of masses larger than 4) or other effects. 

Stellar population analyses do not even agree with these merger rates. Mergers 
or captures of smaller galaxies can occasionally occur, but hierarchical-scheme sub- 
units fusing together and made of gas and stars is not the dominant one by which 
massive elliptical galaxies are made, at least for z < 2 Most massive elliptical 
galaxies hav e a passive evolution since their creation according to stellar popula- 
tion analyses 2 ^ 1 . Mergers are beautiful, spectacular events, but not the dominant 
mechanism by which elliptical galaxies are assembled. Most early-type galaxies with 
a velocity dispersion exceeding 200 km/s formed more than 90% of their current 
stellar mass at redshift z > 2.5 ^21. Elliptical galaxies formed in a process similar 
to monolithic collapse, even though their structural and dynamical properties are 
compatible with a small number of dry mergers 2 ^ 1 , far from the number of mergers 
necessary in our case. Dry mergers do not decrease the galaxy stellar-mass surface 
density enough to explain the observed size evolutional, and the high density of the 
high-z elliptical galaxies does not allow them to evolve into present-day elliptical 
galaxies^. 

Mergers are searched for in the Local Group galaxies too. In the case of the 

Milky Way, some authors try to find evidence of major fusion events of bi g ga laxies, 

but up to now we do not see evidence in favour but against such scenarios^. There 

are minor mergers, of course, and absorption of small clouds of the intergalactic 

medium, but the presence of intermediate mass galaxies at short distances from 

the center, at present or in the past, have yet to be identified. There are certain 

attempts to find something, for instance the recent discovery of a galaxy with a 

relatively large diameter at only 13 kpc from the Galactic center called Canis Major, 

f86l 

but that discovery resulted in a fiascd 2 ^. In any case, two to four major mergers on 

average per galaxy is too much. 

Another argument against the merger scenario of the hierarchical CDM cos- 

1071 

mology is that galaxy formation is controlled by a single parameter 2 -*- 1 . One would 
expect in the merger scenario that the properties of individual galaxies be deter- 
mined by a number of factors related to the star formation history, merger history 
(masses, spins and gas content of the individual merging galaxies), etc., but that 
is not so; all the different parameters of the galaxies are correlated and there is 
only one single independent parameter based on their mass. 

On the other hand, major mergers of disk galaxies of comparable mass should 
give place to elliptical galaxies, so it is not easy to understand in this scenario how 
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the radius of disk galaxies grows. Conselice et alP^l show in fact that there is little 
to no evolution for disk galaxies at z < 1.2, for the K-band, in the stellar-mass 
Tully-Fisher relation, and in the ratio of stellar/total mass. Ferguson et «z|S2] a j so 
find that accretion flows play only a minor role in determining the evolution of the 
disk scalelength. In models in which the main infall phase precedes the onset of star 
formation and viscous evolution, they find the exponential scalelength to be rather 
invariant with time. On the other hand, models in which star formation/ viscous 
evolution and infall occur concurrently result in a smoothly increasing scalelength 
with time, reflecting the mean angular momentum of material which has fallen in 
at any given epoch. 

Furthermore, selection effects go apparently in the opposite direction of observ- 
ing pre-merger galaxies at high redshift. At very high redshift, we are obser ving 
only galaxies with stellar masses over 10 11 M and some of them over 10 12 M ^4 
And at low redshifts there are galaxies of all masses but the average stellar mass is 
much lower than that. Thinking that very massive galaxies at high z are building 
blocks of even more massive low z galaxies is counterintuitive. After 2-4 mergers 
of equal mass galaxies in average, the galaxies should be 4-16 times more massive 
than the original building blocks at high redshift. We should be observing some 
galaxies at low redshift with stellar masses of ~ 10 13 M Q . With a mass-luminosity 
ratio of M*/Ly = 6 (for a very old population ESI if it were younger than 12 Gyr it 
would be lower so the luminosity would be even higher), this would mean galaxies 
with Ly = 2 x 10 12 L 0! y, or absolute magnitude My = —26. Even cD galaxies in 
the centers of the clusters are not as bright as that. Where are these galaxies, then? 

Thinking that low- to intermediate-mass galaxies are the final stage of major 
merger processes is a reasonable possibility, since we cannot see their building blocks 
at high redshift (they are very faint). However, for only high luminous galaxies I 
get more or less the same shrinking factors at high-z with respect to low-z (see 
Fig. [4j. It is not a question of some merging which affects low luminosity galaxies 
more. I could even concentrate our analysis on galaxies with Ly > 1.02 x 10 11 L 0) y 
and, in spite of the poorer statistic, a very strong size e volution between galaxies 
at z = 0.5 — 1 and z > 2 can be appreciated. Refs. I 47 | 55 | even g e ^ higher evolution 



for higher stellar masses, but this has been criticized^. 

Another element that is not consistent with these hierarchical merging scenarios 
is that superdense massive galaxies should be common in the early universe (z > 
1.5), and a non-negligible fraction (1-10%) of them should have survived since that 
epoch without any merging process retaining their compactness and presenting old 
stellar populations in the present universe. However, Trujillo et alP^ find only a 
tiny fraction of galaxies (~ 0.03%) of these superdense massive galaxies in the local 
Universe (z < 0.2) and they are relatively young (~ 2 Gyr) and metal- rich ([Z/H] ~ 
0.2). Clearly a case of how some authors (Trujillo et alP^) try to find proofs in favour 
of the hierarchical merging scenario, and when they find that the observations point 
out exactly the opposite thing of what is expected, instead of claiming that the 
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hierarchical merging scenario is wrong, they try to deviate attention by giving less 
importance to the observed facts for the validation of the standard theory. 

4.4. Quasar feedback 

Fan et al^^ realized that the evolution+merger luminosity solution is not enough 
to explain the strong size evolution and they claimed, "no convincing mechanism 
able to account for such size evolution has been proposed so far" . Nonetheless, they 
have proposed a new mechanism to explain this extraordinary size evolution. Fan 
et ai!^ propose that in elliptical galaxies it is directly related to a quasar feedback: 
part of the energy released by the QSO would be spent to produce outflows of huge 
amounts of cold gas expelled from the central regions, a rapid (few tens of Myr) 
mass loss which induces an expansion of the stellar distribution. 

Although this mechanism might explain some small part of the expansion of 
elliptical galaxies, there are some aspects in this hypothesis which do not fit with 
the observed facts well, at least while we do not clarify some points. First, we do not 
see such supermassive outflows which are necessary to maintain the Fan et al^^ 
idea, although since their life is very short only a few high redshift QSOs would 
show it. Second, elliptical galaxies expend their gas to produce stellar formation 
which gives rise to their stellar mass. If a QSO swept the gas away, no young stellar 
populations would be observed, but there is now compelling evidence for a significant 
post-starburst population in many luminous AGN^^l There is also detection of large 
amounts of warm, extended molecular gas indicating that QSOs have vigorous star 
formation^!, and that the gas is not being expelled. The blue color of host galaxies, 
(B — V) res t ~ indicates a young population too. Third, if QSOs produced 
such massive outflows ejecting the gas of the galaxies, this would also apply to disk 
galaxies. Around 40% of host galaxies in QSOs are disk galaxies^, and flis clear 
that disk galaxies still have gas and active star formation in their disks. Fourth, a 
continuous increase in the average size of ellipticals, as shown in Fig. |3)/left, would 
require the continuous expulsion of gas, but most massive elliptical galaxies have 
had a passive evolution since their creation at z > 2 according to stellar population 
analyses^ 1 , which indicate that the gas was already drained in them at z > 2. Fifth, 
for an average increase of a factor 3 in size in the elliptical galaxies, we would need 
to assume that most elliptical galaxies have hosted very luminous QSOs during their 
lives. With a minimum QSO life tim e of around 40 Myr, as required for the massive 
outflow mechanism of Fan et aZp^, the number of very luminous QSOs should be 
around 1/300 of the number of elliptical galaxies. This number is too large, given 
that the density of QSOs with Lf, i >~ 10 48 erg/s (the necessary luminosity, 5% of 
which is spent to produce outflows larger than 1000 M /yr, as posited by Fan et 
a lM l ~ with M > 2 X 10 10 M Q and R = ±R S dss) is ~ 10" 8 Mpc" 3 ^. 

Nonetheless, I would not dare to say that Fan et al. 's^J hypothesis is incorrect. 
I think it is an elegant and interesting solution to the problem, and it must be 
considered as a serious possibility, provided it is able to solve their caveats, and 
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give some support to the scenario with some observations directly interpretable as 
massive outflows. 



4.5. Rotation or dispersion velocity analysis 

If this hypothesis of lower radius at high redshift for a given mass were true, due 
either to mergers or quasar feedback, we would expect a significant increase in the 
rotation speed or dispersion velocity of those galaxies at high redshift with respect 
to the local ones. For a constant mass within a radius R, one would roughly expect 
rotation speed in a galaxy v ro t oc i? -1 / 2 , and something similar for the dispersion 
velocities in ellipticals. Let us analyze whether this change of velocity is taking 
place. 

I9TI 

An analysis of Marinoni et a/.'s 1 ^ 1 data with galaxies with z < 1.2 does not 
show (see Fig. |6|) a significant change in the rotation velocity-size relationship or 
the r ota tion velocity-absolute magnitude relationship. According to the Saintonge 
et alW^ analysis for low redshift galaxies, the rotation velocity is proportional to 
i?°- 98±0 01 , and is also related to the absolute magnitude in the I-band by their 
eq. (3), so we could derive an average expected velocity from the luminosity of the 
galaxy, v(Mi). In Fig.Hl I see that these relationships remain more or less constant: 
the variations in v ro t/R ' 98 and v rot /v(Mi) are compatible within 1-cr to be null. 
Particularly, the best linear fits give: 

(26.4±5.9) + (8.9±8.8)z, (31) 



#0.98 
Vrot 



= (1.04 ± 0.20)- (0.09 ±0.16>. (32) 



v(Mi) 

Within the error bars, I cannot exclude an increase in these ratios with redshift 
compatible with the hypothesis of radius decrease. An interesting test would be to 
measure the rotational velocity in some of the very compact galaxies with redshift 
3. Although getting a spectrum of these faint galaxies is technically difficult, this 
would give a proof of whether either they are really so compact or the cosmological 
parameters are wrong. 

The critical assumption of a variable effective radius is also counter- argued by the 
proofs in favour of a constant radius^- 1 showing that high redshift first-rank elliptical 
galaxies, with similar absolute magnitudes, have the same vel ocity dispersions as 
low redshift first-rank elliptical galaxies. Van Dokkum et aZ. ^3 point out, however, 
that these galaxies are predicted by galaxy formation models to be those whose 
formation finished at very high redshift and so it would not be surprising that these 
galaxies had the same radius as at low redshift, specifically because they are in 
clusters where mergers may not be likely. 



A surprising new result also points in this directional the velocity dispers ion 
of giant elliptical galaxies with average redshift z w 1.7 from Cimatti et o/.'s^21 
sample is similar or very slightly larger than the dispersion for the same kind of 
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Fig. 6. 

Plot of v ro t/R 98 and v ro t/v(Mi) for the 39 galaxies of the sample from Marinoni 
et a/PH. Squares with error bars are the average of the data (asterisks) with steps 
of 5z = 0.2. The solid line is the best linear fit, given by eqs. (|3"Tj) and ([321) . 

galaxies with the same stellar mass in the local Universe at 240 km/s, while at z = 
it is around 180 km/s. Since Cimatti et a/!^ galaxies at z w 1.7 are more compact 
than the average size with that luminosity at that redshift, it is normal to have a 
slightly higher velocity dispersion. It should be much highe r (o ver 400 Km/s) at 
high redshift since a much lower radius is attributed to them^J, but it is not. This 
result points directly to the conclusion that the galaxies have not strongly changed 
their radii. Other alternative ad hoc ideas in terms of a conspiracy of effects in which 
the dark matter ratio has increased the amount necessary to compensate for the 
radius increase^! sound like a queer coincidence, and have no clear basis in terms 
of galaxy formation scenarios. 

4.6. Discussion on expansion+ evolution models 

All these considerations may make us think that the concordance model cannot 
explain the present data. The other expanding models present similar problems: a 
factor in average size evolution up to z = 3.2 equal to 5.8 for Einstein-de Sitter, 5.6 
for Friedmann with Vt m = 0.3, and 5.9 for the QSSC model. Phenomenologically, it 
is possible to fit the data to any of the expanding models with appropriate evolution 
of galaxies. In practice, this evolution for a constant luminosity galaxy should be 
so strong (up to a factor 200 in average in the luminosity density up to z 3.2; 
systematic errors may change up to a factor 2 this number, but this does not change 
the situation) that the explanations for it seem unrealistic. 

The situation for the expanding models becomes even more dramatic if we go to 
higher redshift. At redshift 6, the linear size of the galaxies assuming a concordance 
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model is even lower, approximately a factor two lower than at z = 3.2 (SB Kg. 4). 
If we were going to consider a reduction in size by a factor 12, all the arguments 
given in this section would become even stronger. 

5. Analysis of the static universe cases 

In the first two static Universe cases, there is an excess of size (~ 20 — 30%) for 
most redshifts. A possible interpretation of this discrepancy is that there may be 
a systematic error in the calculation of the ratio between measured and predicted 
average 6*. In i)3.3l I have discussed the possible systematic errors and I concluded 
that they should be lower than ~ 30 — 40%. These systematic errors might be 
enough to justify the departures of the data with respect to the prediction in Fig. 
[2] for these static models. 

The third case of plasma redshift is much more discrepant and can only fit 
the data with a significant evolution of galaxies, although less strong than in the 
expanding models: a factor 3 instead of a factor 6 in the concordance model for a 
given luminosity from z = 0toz = 3.2. 

5.1. Including extinction 

Apart from the systematic error considerations, another solution to make the two 
first static cosmological models compatible with our data would be related to ex- 
tinction rather than the evolution. Extinction would make the galaxies look fainter, 
which means that, through Eqs. ([2]) and (j4j, when the corrections of extinctions are 
made, the luminosity is larger and (9» is smaller than their values without correc- 
tions. Let us check this hypothesis with a rough calculation. Instead of Eq. ([5]), the 
inclusion of the IGM extinction with absorption coefficient K (area per unit mass) 
will give the following relationship 

/'d J] , r , 1 + Z (d A ) i 

^ ° irs!v wl (33) 

I have assumed a constant dust density pdust along the line of sight, which is an ap- 
propriate approximation for a homogeneous Universe without expansion and moder- 
ate amounts of dust ejection by the galaxy. If we considered an expanding Universe 
with a strong dust emission rate by the galaxies, we should include pdust(r) within 
the integral, but it is not the case here. The absorption coefficient can approximately 
be described with a wavelength dependence: 

k(\) = k(\ v )(J- > ) a . (34) 

Hence, and using Eqs. (fT5j) and (fi~7|) , 

L v , rest = AnF Krest d 2 L e^^ [{1+z)m - {1+zra \ (35) 
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with tyl — 1 for the cosmology with linear Hubble law, and m — for the tired light 
case, ay = n(\y) Pdust is the absorption in V per unit length, which means there 
are 1.086ay magnitudes in V of extinction per unit length. 

The value of the exponent a is not totally independent of A but I take it ap- 
proximately as constant. I adopt a = 2, as observed in near-infrared bands in our 
Galaxy^Sl. For lower wavelengths (optical, ultraviolet) a would be lower. Since most 
of the sources have the wavelength equivalent to V-rest in the near-infrared, and 
the extinction c urve of dust in the intervening QSO absorbers resembles the SMC 
extinction curve^^, this approximation is reasonable. 

The absolute value of ay is not well known and neither do we know whether 
it is significant or null. There are only some constraints for the maximum value 
(e.g-J ^-H although based on standard cosmology). The values of ay which give the 
best fit to our data are: ay — 1.6 x 10~ 4 Mpc" 1 for the linear Hubble law case, 
and ay = 3 .4 x 10 -4 Mpc -1 for the simple tired light case. Assuming n(\y) ~ 10 5 
cm 2 /gr^2l, the value for the dust density necessary to produce such an extinction 
would be pdust ~ 6x 10~ 34 g/cm 3 , and pdust ~ 1.2 x 10 -33 g/cm 3 respectively, 
which is within the range of possible values. Vishwakarm gi ves va lues of pdust = 
3 — 5 x 10~ 34 g/cm 3 , but for the QSSC model. Inoue & Kamaya^^ allow values as 
high as pdust ~ 10 33 g/cm 3 for the high z IGM within the standard concordance 
cosmology. F or co mparison, the average baryonic density of the Universe is (taking 
il b = 0.042; EEH) Pb = 3.9 x 10~ 31 g/cm 3 , so this would mean that IGM dust 
constitutes 0.15 or 0.30% of the total baryonic matter, reasonable amounts. 

Whether the extinction used in the models with extinction would be grey or 
would introduce some small reddening is not totally clear, but this is not a question 
for the present analysis. I just note that reddening in the optical would depend on 
the value of a in the visible at intermediate to high redshift and the variability of 
k(X) with respect to A in the UV. The features of dust extinction in the UV are not 
easy to model in an IGM with unknown composition. 

With this simple correction for extinction, the results are significantly improved, 
as shown in Fig. [7] The tired light case gives a better fit. 



5.2. Comparison with angular size test for ultra-compact radio 
sources 

Compact radio sources have been used by several authors to carry out the angular 
size test because these sources were thought to be free of evolutionary effects. How- 
ever, the different results obtained with these sources has raised the suspicion that 
they may not be such good standard rods. Apparently, these rods are somewhat 
flexible. For example, KellermanrPSl claimed that the angular size test for these 
sources fitted Einstein-de Sitter expectations ver y well, when Einstein-de Sitter 
was the fashionable model. Jackson & Dodgsoiff^ claimed the opposite: that it was 
not compatible with Einstein-de Sitter, and that, given that VL m = 0.2, the best 
fit for the cosmological constant term was Cl\ — —3.0; flat cosmological models 
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Fig. 7. Log— log plot of the average of log 10 #*, where 0* is the equivalent angular size, vs. redshift 
(z). Bins of Alog 10 (^) = 0.10. The two plots are for the first two static cosmologies with the 
inclusion of a constant IGM extinction which gives the best fit. 



were excluded with > 70% C.L. JacksorffSI, in the era of the concordance model 
as the fashionable cosmology, again carried out the analysis of the same data used 
by Jackson & Dodgso doing some new corrections due to selection effects and 
bias, and they get the best fit for f2 TO = 0.29, J7a = 0.37, compatible within 1-a 
with the concordance model. With further data, Jackson & JannettsP^ get the best 



fit for fL 



0.25 



-0.04 
-0.03' 



n A = 0.97 



+0.09 
-0.13 



i C.L.). It seems that the general trend 



is to obtain the result expected from fashionable cosmologies on the date in which 
the test is carried out, and when incompatibilities appear, some selections effects, 
biases, small evolution effects are sought to try to make the results compatible. In 
my opinion, this is not a very objective way to do science, but let us leave the 
discussion of the methodology of cosmology aside. 

One important selection effect is derived from the fact that linear sizes depend 
on radio luminosities. JacksorP^ tries to take this effect into account and suggests 
a method of correcting it: binning the data into groups of 42 points in the redshift 
distribution and taking as representative of each bin the mean of points between 
11 to 17 within each one, counting from the smallest objects. This is supposed to 

— 1/3 

be done to compensate for the dependence R ra d oc £ rad for a given redshift and 
for the fact that the lowest luminosity points cannot be observed at high redshift. 
In my opinion, this is not the right way to correct the selection effect. JacksorP^l is 
just doing a median which gives more weight to smaller objects, but this median is 
shifted at high redshift by the lack of low luminosity objects. 

Nevertheless, my concern is not about Jackson's method of correcting for the 
Malmquist bias, but about the relationship between radius and luminosity. The 

II Ql 

relationship in Fig. 1 of Jackson is applicable only to the concordance model and 
it will be far different for static Universes. Particularly in Figs. 1 and 2 of JacksorP^ 
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we see that the linear size is almost independent of redshift between z = 0.5 and 
z = 4: around 10 pc, with some scattering due in part to the range of luminosities for 
each redshift. However, if I use eq. 1151) of the static Universe with linear Hubble law 
instead of eq. Q of the concordance model to calculate the angular size distances, I 
see that the linear sizes at z « 3.5 should be a factor of 10 larger (~ 100 pc instead of 
~ 10 pc). Therefore, the linear sizes would not be independent of redshift but highly 
dependent on it. In such a case, Figs. 1 and 2 of JacksoiffSI would be transformed 
into a plot with a continuous increase in linear size with radio luminosity for all 

1/2 

ranges of luminosity for all rcdshifts: roughly R ra d oc L d . Since at high redshift we 
cannot observe objects with low L ra d, this means that we would be losing objects of 
low Rrad at high z, so the median of 9 would be highly overestimated with respect 
to low to intermediate redshift objects. This would explain why the angular size 
test for ultra-compact radio sources does not give a z _1 dependence. Therefore, 
a static Universe is not excluded by this test unless we demonstrate with data at 
low redshift that the radii of the ultra-compact sources does not depend on their 
luminosities. 

Ultra-compact objects could be used to carry out a right angular size test, but 
we cannot directly compare objects of low luminosity at low redshift with objects 
of high luminosity at high redshift because we have no guarantee that the linear 
size is the same in both cases (and it is not valid to assume a cosmological model 
a priori to prove that it is good a standard rod, because the method should be 
independent of any cosmological assumption if we want to derive from it which 
is the best cosmological model). We should either i) compare objects of the same 
luminosity (different for each cosmology), or ii) define a as in the present paper 
in which we need to calibrate the size-luminosity relationship in the low-z Universe. 
This second option has the caveat that we do not have very high luminosity compact 
radio sources at low redshifts so we need to extrapolate the local radius-luminosity 
relationship for high luminosities. 

5.3. Hubble diagram for the different cosmologies 

In Fig. [SI I show the different distance moduli for the different cosmologies without 
extinction, together with some real data of SNe la compiled by Kowalski et al\ ^^ \ 

m Vtres t - M Vtres t = 5 log 10 [d L (z) (Mpc)] + 25. (36) 

In the first two static models, if we wanted to include the extinction, we should sum 
A v , rest (z) = g°^ c + ^ [(1 + z) m - (1 + z)- a ] to the distance modulus. 

One aspect is remarkable: the value of the distance modulus for the concordance 
model is very similar to its value for the static model with a linear Hubble law, and 
it can be seen in Fig. [8J how the data of SNe la are approximately compatible with 
this scenario. The fit for the concordance model over the data gives a reduced x 2 : 
Xr = 3.34; while for the static model with linear Hubble law \r = 4.20, slightly 
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worst but not by much: the concordance model reduces the % 2 by only 20%. LerneJ^Sl 
also show this by comparing the residuals of the Hubble diagram in the concordance 
model and in the static universe ansatz and realizing that they are both similar. Is it 
a coincidencdB? With a slight extinction of ay ~ 1 — 2 x 10~ 4 Mpc~ 3 (the range of the 
best fit obtained in i)5.1[) . the agreement is also quite conspicuous at least for z < 1. 
For the simple tired light model, only with extinction of the order ay ~ 5 x 10 
Mpc~ 3 (not very far from the best fit obtained in fc|5.1[) would get the coincidence for 
z < 1. For the plasma redshift tired light model without extinction, the agreement 
with SNe la for z < 1 is also acceptable, as noted by BrynjolfssorpSl. This means 
that, with the static models, we can fit nearly the same Hubble diagrams as the 
concordance model with its cosmological constant, particularly for supernovae fits. 
It is not the place here to extend the discussion on the analysis of the compatibility 
of the static model with a linear Hubble law and the supernovae la data; this would 
require a discussion of the systematic errors, the selection effects, etc. At present, I 
just want to emphasize that there are no major problems to make compatible the 
static model of linear Hubble law with SNe la data. 

For other cosmologies, and without extinction, the difference from the con- 
cordance model is larger, especially for the highest redshifts. This does not 
mean that they are discarded, because the objects used as standard can- 
dles (particularly supernovae) might have an absolute magnitude which is not 
consta nt with redshift or some extincti on along their lines of sight. Several 
authora^U- 06 ! 107 ! 108 ! 109 ! 110 ! 111 ! 32 ! ^ also think that SNe la data used to de- 
rive J7a = 0.7 are affected by several systematic uncertainties that make the A-CDM 
cosmology uncertain. 



5.4. Evolution, and other tests 

Another topic to discuss is the apparent evolution of galaxies at different redshift. 
I have not included any evolution correction for the static models in the plots of 
Fig. [7] although some slight evolution might be compatible with them, in the sense 
that brighter populations are present at higher redshifts. Some evolution might be 
present because, as discussed in ^3. 2\ elliptical galaxies get lower angular sizes at 
high redshifts than disk galaxies, possibly due to the different mean ages of their 
populations or merger rate; although some of these differences might be due to sys- 
tematic errors too, as said in M3.31 Also, the most massive galaxies present a higher 



b The degree of coincidence depends on the maximum redshift, z max , we use. For instance, the 
value of f?A in a flat Universe which gives a best minimum square fit to the Hubble diagram for a 
static model with linear Hubble law and without extinction is Ha ~ 0.39+0.108z m aa; — 0.00852^ x . 
For Zmax = 6, as in our plot of Fig. [8] it gives £7a = 0.74, but it falls to 0.55 if we only consider 
it up to redshift 1.7, as usually for SNe la data. If we set f2 m = 0.3 and we searched for the best 
value of f2j\ with any value of Qtotali we would get the best fit for fs 0.76 + 0.171z ma:E — 
0.053x^2. + 0.0038,2^2.. For z max = 6, it gives the best fit for Qa = 0.70, and it increases to 0.91 
if we only consider it up to redshift 1.7. 
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Fig. 8. 

Distance modulus as a function of redshift for different cosmological models. Data 
of supernovae la compiled by Kowalski et aZP^l wer e added. Right plots are zooms 
of the left plots. 



ratio of angular size with respect to local which may be interpreted 

here as a real higher compactness with respect to the least massive galaxies. Note, 
however, that, as said above, this extra-compactness can be understood in terms 
of fluctuations due to noise preve ntin g the recovery of the extended low surface 
brightness halos in the light profild^. 

There is also a slight color at rest evolution^ and a mass-luminosity ratio 
evolution, as said in £14.21 but these mass calculations are subject to important 
errors depending on the synthesis model, IMF assumption, etc.; hence, there is 
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a wide range of possible values of mass-luminosity ratio evolution. The bluer (at 
rest) color of high redshift galaxies might be due to bias, because bluer galaxies with 
younger populations are brighter. We must also bear in mind that photometric errors 
are larger at higher redshift (because of the fainter fluxes), the photometric redshifts 
have higher uncertainties and consequently might affect the deter minati on of the 
luminosities at rest, etc. At present, with the data used from Refs. | 24 | 41 | a t z > i 5 
the correlation between z and {B — V) rest is: (z(B — V) rest ) — (z)((B — V) rest ) = 
-0.048 ± 0.021 for elliptical galaxies and -0.033 ± 0.015 for disk galaxies. The 
variation in color with the variation in redshift is much lower than the dispersion of 
colors and possible systematic effects. Therefore, we cannot exclude the possibility 
that the luminosity evolution is small enough to be within the error bars of our 
data with the static models. 

Note however that the evolution of some quantities per unit comoving volume 
[for instance, the star formation ratio, expressed in mass per unit time per unit 
comoving volume, is claimed to be significantly higher in the paslP^j must be 
corrected with respect to the concordance cosmology with Eqs. (f22|) . (|24[) . which 
reduces by a factor of 10 at z — 2 the ratio in the static model with a linear Hubble 
law with respect to the concordance one. 

In any case, there may be a real evolution which should be explained, either 
in expanding or static models. Note that "static" does not mean "no evolution"; 
"static" means "no expansion" ; there is not necessarily a contradiction in observ- 
ing evolution in a static Universe, provided that the creation of galaxies is not a 
continuous process. 

Explaining the cosmic microwave background radiation (CMBR) and its 
anisotropics is not the purpose of the present paper. Note, however, that there exist 
alternative scenari os to explain them apart from the model of s tandard hot Big 
Bang (e.g., I * \ see the discussion in another pape r of min §1). Concerning 
other arguments/tests in favor of the expansion (e.g., 14411161117^ 

, we must bear in 

mind that they are usually a matter of discussion. Tests such as the time dilation in 
SNe la, which were claimed to be a definitive proof of t he expansion of the Universe, 
find counterarguments and criticisms from who claim that 

a static Universe is compatible with the data. The same thing happens with any 
other test, including the present one of angular size. Apart from the present angular 
size test, there are other tests that also p resent results in fayor of a static Universe 
and against an expanding Universe (e.g., I 12 ! 121 ! 122 ! 123 ! 124 ! 125 ! 32 ! 45 ' 46 !). 

Perhaps the most immediate problem with the static Universe is understanding 
the cause of the redshift of galaxies, but there are several proposals for alternative 
mechanism to produce redshifts without expansion or Doppler effect, so the hypoth- 
esis of a static Universe is not an impossible one. Other facts, such as the formation 
of the large scale structure, the creation of the light elements, etc., also provide 
alternative explanations different from the standard model. Fu rther discussion of 
all the questions raised in this paragraph are given in my review^SI. 
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5.4.1. UV surface brightness test 

The main discrepancy in the different tests such as angular size, surface brightness, 
Hubble diagram, etc., is the evolution of galaxy luminosities, which is very large 
for the defenders of the expansion and not so large for the defenders of the static 
models. LerneJ3Sl proposes a test of the evolution hypothesis that is also useful in 
the present case. There is a limit on the UV surface brightness of a galaxy, because 
when the surface density of hot bright stars and thus supernovae increases large 
amounts of dust are produced to absorb all the UV except that from a thin layer. 
Further increase in surface density of hot bright stars beyond a given point just 
produces more dust, and a thinner surface layer, not an increase in UV surface 
brightness. Based on this principle, there should be a maxim um s urface brightness 
in UV-rest wavelengths independent of redshift. Scarpa et a measured in low 
redshift galaxies a maximum FUV (1550 A at rest) emission of 18.5 magAs/arcsec 2 
(the average is 24 — 25 magAs/arcsec 2 ); no galaxy should be brighter per unit 
angular area than that. 

Using eqs. (|2"T1) . (|2"51) with the flux at wavelength 1550 A at res1@ from the 
subsample MS 1054-03 in Trujillo et al¥^ galaxies, and the angular sizes @ 155Q ^ = 

1.146V ^ for n s > 2.5 and 6 ^ = 1.106V ^ for n s < 2.5, I get the values 
of SBq for all the galaxies in this subsample for the expanding or static cases 
(Fig. [9]). For the expanding universe, many galaxies have average intrinsic surface 
brightness (SBq) lower(brighter) than 18.5 magAs/arcsec 2 , the galaxy MS 1054- 
03/1356 being the brightest one per unit angular area: 14.8 magAs/arcsec 2 (30 
times brighter than the limit). The angular size of this galaxy MS 1054-03/1356 
(0.027" circularized in V; 0.031" in FUV) might be affected by some error since it 
is below 0.125" in V, but even an error of 100% in angular size would produce an 
error of 1.5 magnitudes in surface brightness, not 3.7 magnitudes as we observe here. 

The dispersion (r.m.s.) of ^ 50A might be around 20% Fig. 2), which means an 
uncertainty of around 0.4 mag/arcsec 2 (1-c) in SBq, much lower than the differences 
between SBq and its limit of 18.5. Moreover, even avoiding galaxies with angular 
size less than 0.125", there are some galaxies with surface brightness over the limit, 
up to 6 times brighter than the limit. Too high intrinsic FUV surface brightness. 
LerneiS^l also argues why other alternative explanations (lower production of dust 
at high redshift, winds or others) are not consistent. However, for the static models, 
all galaxies have average intrinsic surface brightness (SBq ) within 9 higher (fainter) 
than the limit 18.5 mag^s/arcsec 2 , a result which may be interpreted again in favor 
of the static scenario. 



c An interpolation with the publicly available fluxes in filters U, B, V, Isi4 and J was used to get 
it. These fluxes are given in flux per unit frequency so the dimming factor must be multiplied by 
a factor (1 + z) with respect to the flux in the whole filter. That is, SBq = SB(1 + z) 3 for the 
expanding case, and SBq = SB for the static case. 
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Fig. 9. Intrinsic FUV-rest surface brightness of the galaxies in the subsample MS 1054-03 in an 
expanding (left) or static (right) universe. The dashed line stands for the minimum value of 18.5 
mag/arcsec 2 over which all galaxies should be located. Points with circles stand for data with 
9 < 0.125". 



5.5. Is a static model theoretically impossible? 

Apart from the discussion on the observations, which are inconclusive, a static 
model is usually rejected by most cosmologists on the grounds of a belief/prejudice 
that a static model is impossible. However, from a purely theoretical point of view, 
without taking into account the astronomical observations, the representation of the 
Cosmos as Euclidean and static is not excluded. Both expanding and static space 
are possible for the description of the Universe, even with evolution. 

Before Einstein and the rise of Riemannian and other non-Euclidean geometries 
to the stage of physics, attempts to describe the known Universe with a Euclidean 
Universe were given, but with the p roblem of justifying a stable equilibrium. Within 
a relativistic context, Einstei n 1 proposed a static model including a cosmologi- 
cal constant, his biggest blunder according to himself, to avoid a collapse. This 
model still has prob lems to guarantee the stability, but it might be solved somehow. 
Narlikar & Arp^S so lve it within some variation of the Hoyle-Narlikar conformal 
theory of gravity, in which small perturbations of the flat Minkowski spacetime 
would lead to smal l oscillations about the line element rather than to a collapse. 
Boehmer et al\ ^^ analyze the stability of the Einstein static universe by consid- 
ering homogeneous scalar perturbations in the context of f(R) modified theories 
of gravity and it is found that a stable Einstein cosmos with a positive cosmolog- 
ical const ant is po ssible. Other authors solve it with the variation of fundamental 
constants I 133 | 134 | Another idea by Van FlanderrP^S] j s that hypothetical gravitons 
responsible for the gravitational interaction have a finite cross-sectional area, so that 
they can only travel a finite distance, however great, before colliding with another 
graviton. So the range of the force of gravity would necessarily be limited in this 
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way and collapse is avoided. ^^^^ 
As said in £14.1[ the very concept of space expansion has its own problemJS^IMl 
The curved geometry (general relativity and its modifications) has no conservation 
of the energy-momentum of the gravity field (the well-known problem of the pseudo- 
tensor character of the energy- momentum of the gravity field in general relativity). 
However, Minkowski space follows the conservation of energy-momentum of the 
gravitational field. One approa ch w ith a material tensor field in Minkowski space is 
given in Feynman's gravitatio n^*^ , where the space is static but matter and fields 
can be expanding in a static space. It is also worth mentioning a model related 
to modern relativistic and quantum field theories of basic fundamental interactions 
(strong, weak, electro- magnetic): the relativ istic field gravity theory and fractal 

[Tool 

matter distribution in static Minkowski spac e 1 

Olber's paradox for an infinite Universe also needs subtle solutions, but ex- 
tinction, absorption and reemission of light, fractal distribution of density and the 
mechanism which itself produces the redshift of the galaxies might have something 
to do with its solution. These are old questions discussed in many classical books 
on cosmology (e.g.,222, ch. 3) and do not warrant further discussion here. 



6. Conclusions 

Summing up, the main conclusions of this paper are the following: 

• The average angular size of the galaxies for a given luminosity with redshifts 
between z = 0.2 and 3.2 is approximately proportional to z~ a , with a between 
0.7 and 1.2, depending on the assumed cosmology. 

• Any model of an expanding Universe without evolution is totally unable to fit 
the angular size vs. z dependence. The hypothesis that galaxies which formed 
earlier have much higher densities does not work because it is not observed here 
that the smaller galaxies are precisely those which formed earlier; in any case, 
the galaxies observed today were formed mostly at redshifts not very different 
from the galaxies observed at higher redshifts. A very strong evolution in size 
would be able to get an agreement with the data but there appear caveats 
to justify it in terms of age variation of the population and/or mergers and/or 
ejection of massive outflows in the quasar feedback. An average of the necessary 
two to four major mergers per galaxy during its lifetime is excessive, and neither 
is it understood how massive elliptical galaxies may present passive evolution in 
this scenario or how spiral galaxies can become larger during their lifetimes. The 
depletion of gas in ellipticals by a QSO feedback mechanism does not appear to 
be in agreement with the observed star formation and other facts. Moreover, no 
evolution is observed in the rotation/dispersion velocities and, the FUV surface 
brightness turns out to be prohibitively high in some galaxies at high redshift. 

• Static Euclidean models with a linear Hubble law or simple tired light fit the 
shape of the angular size vs. z dependence very well: there is a difference in 
amplitude of 20-30%, which is within the possible systematic errors. An extra 



2, 2010 16:14 WSPC/INSTRUCTION FILE sizes 



Angular size test 37 

small intergalactic extinction may also explain this difference of 20-30% . Some 
weak evolution of very high redshift sources is allowed, although non-evolution 
is a possible solution too. For the plasma redshift tired light static model, a 
strong (albeit weaker than in expanding models) evolution in galaxy size is 
necessary to fit the data. The SNe la Hubble diagram can also be explained in 
terms of these models. 
• It is also remarkable that the explanation of test results with an expanding 
Universe requires four coincidences: 

(1) The combination of expansion and (very strong) size evolution gives nearly 
the same result as a static Euclidean universe with a linear Hubble law 
alone: 6 cx z^ 1 . 

(2) This hypothetical evolution in size for galaxies is the same in normal galax- 
ies as in QSOs, as in radio galaxies, as in first ranked cluster galaxies, as 
the separation among bright galaxies in clusters. Everything evolves in the 
same way to produce approximately a dependence 9 cx 

(3) The concordance model gives approximately the same (differences of less 
than 0.2 mag within z < 4.5) distance modulus in a Hubble diagram as 
the static Euclidean universe with a linear Hubble law. 

(4) The combination of expansion, (very strong) size evolution, and dark mat- 
ter ratio variation gives the same result for the velocity dispersion in ellip- 
tical galaxies (the result is that it is nearly constant with z) as for a simple 
static model with no evolution in size and no dark matter ratio variation. 

These four coincidences might make us think that possibly we should apply 
Occam's razor "Entia non sunt multiplicanda praeter necessitatem", we should 
use the simplest models that can reproduce the same things as a complex model 
with many more free parameters does. 

It would be an irony of fate that, after all the complex solutions pursued by 
cosmologists for the last century, we had to come back to simple scenarios such as 
a Euclidean static Universe without expansion. None the less, we cannot at present 
defend any of these simple models apart from the standard one because this would 
require other analyses. The conclusion of this paper is just that the data on angular 
size vs. redshift present some conflict with the standard model, and that they are 
in accordance with a very simple phenomenological extrapolation of the Hubble 
relation that might ultimately be linked to a static model of the universe. 
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